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Abstract 

We study KKLT type models with moduli-mixing superpotential. 
In several string models, gauge kinetic functions are written as linear 
combinations of two or more moduli fields. Their gluino condensation 
generates moduli-mixing superpotential. We assume one of moduli 
fields is frozen already around the string scale. It is found that Kahler 
modulus can be stabilized at a realistic value without tuning 3-form 
fluxes because of gluino condensation on (non-)magnetized D-brane. 
Furthermore, we do not need to highly tune parameters in order to 
realize a weak gauge coupling and a large hierarchy between the grav- 
itino mass and the Planck scale, when there exists non-perturbative ef- 
fects on D3-brane. SUSY breaking patterns in our models have a rich 
structure. Also, some of our models have cosmologically important 
implications, e.g., on the overshooting problem and the destabiliza- 
tion problem due to finite temperature effects as well as the gravitino 
problem and the moduli problem. 
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1 Introduction 



String/M theory is a promising candidate for unified theory including grav- 
ity. Its 4D effective theory, in general, includes many moduli fields. Their 
vacuum expectation values (VEVs) play an important role in particle physics 
and cosmology. That is because those VEVs determine coupling constants 
such as gauge and Yukawa couplings, and physical scales like the Planck 
scale Mp and the compactification scale. How to stabilize moduli (at re- 
alistic values) is one of important issues to study in string phenomenology 
and cosmology. How to break supersymmetry (SUSY) is another important 
issue. It is expected that nonperturbative effects, which fix moduli VEVs, 
may also break SUSY. A number of studies on these issues have been done, 
and some of those vacua correspond to anti-de Sitter (AdS) vacua. 

Recently, flux compactification is studied intensively, because it can sta- 
bilize some of moduli. In type IIB string models, complex structure moduli 
and the dilaton can be stabilized around the string scale Mstring = a'"^^^, 
which is expected to be of 0(10^^) - 0(10^^) GeV, but Kahler moduh fields 
are not stabilized [1]. On the other hand, in type IIA string models all of 
moduli can be stabilized [2], and in hctcrotic string models complex struc- 
ture and volume moduli can be stabilized, but the dilaton VEV is not fixed. 
However, in the heterotic case, the compact space is generically non-Kahler, 
and it is mathematically difficult to treat it by present knowledge.^ 

Within the framework of flux compactification in type IIB string models, 
a simple model has been proposed in Rcf. [4], where the remaining Kahler 
modulus T is stabilized by nonperturbative effects such as gluino condensa- 
tion. Still, the potential minimum corresponds to the SUSY AdS vacuum 
at the first stage. Then, an anti-D3 (D3) brane is introduced (at the tip of 
throat) in order to uplift the vacuum energy and realize a dc Sitter (dS) or 
Minkowski vacuum. It shifts the potential minimum, and breaks SUSY in a 
controllable way. That is the KKLT scenario. 

Furthermore, soft SUSY breaking terms in the KKLT scenario have been 
studied in Ref. [5, 6], and it has been found that in the KKLT scenario 
the F-term of modulus field T is of 0{m^/2Mp/4:7r'^), where is the 

gravitino mass. That is, the modulus F-term F'^ and the anomaly medi- 
ation [7, 8] are comparable in soft SUSY breaking terms. That leads to a 
quite novel pattern of SUSY breaking terms. (See for their phenomenolog- 
ical aspects Ref. [9, 10, 11].) It is useful to introduce the parameter a as 

''Recently, a solution that compact space is conformally Kahler is found in Ref. [3]. 
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a = /jji ^, ■) [9]) ill order to represent the ratio between the modulus 

mediation and the anomaly mediation. One of phenomenologically inter- 
esting features in the modulus-anomaly mixed mediation is the appearance 
of a mirage messenger scale [9], where the anomaly mediation at the GUT 
scale can cancel the renormalization group effect under a certain condition. 
The mirage messenger scale is estimated as {m^/2/Mp)"^^MGUT- 
That is, soft SUSY breaking terms appear as the pure modulus mediation 
at the mirage scale. The mirage scale depends on the ratio between the 
modulus i^'-term and 171^/2 ■ The original KKLT model leads to a ~ 1, and 
the mirage scale is the intermediate scale. 

Moreover, it has been pointed out in [12] that if the mirage scale is 
around 0{1) TeV, such model has an important implication on solving the 
little SUSY hierarchy problem. At any rate, the mirage scale is determined 
by the ratio a, which has phenomenologically interesting implications. 

The KKLT model is also interesting in cosmology. The anomaly me- 
diation is sizable, that is, the gravitino mass is 0(10) TeV. The modulus 
is much heavier. Thus, we may avoid the cosmological gravitino/moduli 
problem [6, 9, 10, 11]. 

In this paper, we study a modified KKLT scenario. A gauge kinetic 
function is a mixture of two or more moduli fields in several string models, 
e.g., weakly coupled heterotic string models [13, 14], heterotic M models 
[15, 16, 17, 18. 19]. type IIA intersecting D-brane models and type IIB mag- 
netized D-branc models [20, 21]. Suppose that the gauge kinetic function / is 
a linear combination of the dilaton S and the modulus T like / = mS + wT. 
Following Ref. [4], we assume that S is frozen already around M string ^- We 
also assume gluino condensation generates nontrivial T-dependent superpo- 
tential, which form is expected to be e~'^^ ~ Be^^^ . The original KKLT 
model corresponds to the case that B = 0{1) and b is positive. However, in 
our case the constant B can be very suppressed depending on mc{S), and 
also the exponent coefficient b can be negative even though it is generated 
by asymptotically free gauge sector ^. We study moduli stabilization and 
SUSY breaking with such potential terms, by adding uplifting potential. 
One feature of such SUSY breaking is that the ratio between F'^ / {T + T) 

^Alternatively, both of them may remain light in some models. We will study such 
models separately in Ref. [22]. 

®A similar superpotential like e*"^ with 6 > has been studied for non-asymptotically 
free models in Ref. [23]. However, we stress that our superpotential is generated by 
asymptotically free gauge sector. 
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and m3/2 can vary by a value of B. That is phenomenologically interesting, 
because we would have richer structure of SUSY spectra. Moreover, the 
superpotential term with b < has important implications on cosmology. It 
may avoid the overshooting problem and the destabilization problem due to 
finite temperature effects. 

This paper is organized as follows. In section 2, we review on flux com- 
pactification and the KKLT model. In section 3, we give concrete string/M 
models, where gauge kinetic functions depend on two or more moduli fields. 
Then, in section 4 we study moduli stabilization and SUSY breaking. In 
section 5, we discuss implications on SUSY phenomenology and cosmology. 
Section 6 is devoted to conclusion. In Appendix A, detailed analysis on the 
potential minimum is summarized. 



2 Review on flux compactiflcation and KKLT model 
2.1 Flux compactification 

In this subsection, we review on flux compactification and in the next sub- 
section we review on the KKLT model. 

We consider the Type IIB 03/07 oricntifold 4D string model on a 
warped Calabi-Yau (CY) threefold with /ii^i(CY) = 1. In addition, we 
introduce RR and NSNS 3-form fluxes F^^ and H^^, which should be 
quantized on compact 3-cycles C3 and C3 such as 

^ ^ F3«« G 27rZ, f Hi"^ G 27rZ. (1) 



H^^ A Fg^^ + Q^"^"^ = 0, (2) 



27ra' Jcs 27ra' Jc^ 

Furthermore, these fluxes should satisfy the RR tadpole condition 
1 

(2^a'2 Jm, 

where is the RR charge contribution of local objects including D3- 

brane, wrapped D7-brane and 03-planes in the D3-brane charge unit. In 
this flux compactified type IIB string model, we can fix dilaton and complex 
structure moduli including a warp factor but not the Kahler modulus around 
the string scale. Thus, only the Kahlcr modulus remains light. 

We are interested in two moduli, that is, one is the dilaton S and an- 
other is the overall Kahler modulus T, although the dilaton is frozen around 
M string- They are given by 

27rS = e-"^ - ico, 2'kT = v^J.^ - 104. (3) 
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Here, is lOD dilaton and ve is a volume of CY in the string unit within 
the Einstein frame. The lOD Einstein metric is given by the string metric 
9m^n^ = e'^^'^gfiN- Thus, the CY volume in the string frame v is written 
hy V = e^'^f'^VE- The axionic mode cq is the RR scalar and C4 is the 4D 
Poincare dual of (1,1) part of 4-form RR potential. For example, these 
fields are related to gauge kinetic functions on D3 and non-magnetized D7 
brane fop 

fD3 = S, fD7 = T, (4) 

where (Refop) = g]jp and gop is the gauge coupling on the Dp-brane. In 
this perturbative description, it is natural that (ReS), {ReT) = 0(1). In 
addition, these VEVs are related to the physical scales such as the 4D Planck 
scale Mp and the compactification scale Mkk [16, 6] 

where MxR/Mstring = {v^^/'^). We need the condition (ReT) > (ReS). 

With the 3-form flux and orientifold planes in the compact space, the 
dilaton S and the complex structure moduli arc frozen around Mstring and 
a background metric is warped , then the metric in the lOD Einstein frame 
is given by [1] 

dsl^ = -^e^^'^y^g^dx^^dx'^ + e-^''^^'^ v^;^ gmndy^ dy^ . (6) 

gs 

Here, gg = e^'^^ is string coupling, gj^^ is 4D Einstein metric and gmn is 
unwarped compact CY metric which is normalized as / d^y^/detgmn = 
(27ra'^/^)^. The y™ dependence of the warp factor e'^^^y^ on the throat 
is studied in [24, 25] and in generic point we have e^^^^^ ~ 1. A minimum 
of warp factor can be treated as complex structure deformation from CY 
conifold, and the warp factor on the tip of the throat e"^"'*" can be stabilized 
by 3-form flux such that [1, 24], 

A 

ao = exp[(-27r/i)/(35./)] = (7) 
Here, /, h are given by RR and NSNS 3-form fluxes 

^-^ = d^/.^- -2.k=^lHi'^ MeZ, (8, 
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where A is 3-cycle and B is dual cycle of A near conifold singularity. Thus, if 
h ^ gsf , we can produce exponentially large hierarchy in this string model, 
like the Randall-Sundrum model [26]. 

We consider moduli stabilization within the framework of 4D A'^ = 1 
effective supergravity. Here and hereafter, we use the unit that Mp = 
2.4 X 10^^ GeV = 1. Here we neglect the warp factor dependence of po- 
tential, because warping effects are not important in generic point of the 
compact CY space except for the small region on the warped throat. Hence, 
since moduli are bulk fields, they may not be affected by warping. The 
stabilization of S and the complex structure moduli C/" is as follows. The 
3-form flux in the compact space Mq generates the following superpotential, 
[27] 

Wfiux = W{S,U")= f Gs A (9) 

JMe 

where G3 = — 2mSH^^ and Q is the holomorphic 3-form on CY. With 
the following Kahler potential, 

K = -ln(5 + 5) -31n(r + T) -ln(-i / QaTi), (10) 
the scalar potential is written as 

where DaW = {daK)W+daW, = dad^K, a, b are summed over all mod- 
uli fields and i,j are moduli fields excluding T because of no-scale structure. 
The potential minimum, i.e. the F-flatness condition, is obtained as [1, 28] 

DucW = [ G3 A Xa = 0. 

Here, we have used djjaVt = — {du^K) Q.+Xai where Xa is a basis of primitive 
(2,1) forms. ^ Thus, these moduli are generically stabilized at values of order 
unity around M string but the Kahler modulus cannot be stabilized at this 

'^The primitivity means g^'^Xijk = 0- 
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(12) 
(13) 



stage. The stabilization of the dilaton and the complex structure moduli 
is possible, because the degree of freedom of G3 is equal to the number of 
the dilaton S and complex structure moduli, that is 2(1 + /i2,i)- From (12) 
and (13), the 3-form flux G3 can be decomposed to (0,3) form G(o,3) ^^id 
primitive (2,1) forms G^2i)- ^ result, G3 is imaginary self-dual (ISD), 
such that *qGs = iG^. That leads to 

(27r)V^iVfl^ ^ / Hi^'^Fi^^='^ [ GsAGS, 

= ^ / d^'yV^GmniGT'"' >0, (14) 

2 ■ 3! J Me 

then from (2) we can see that we actually need negative charge contribu- 
tions, such as 03-plane, wrapped D7-brane or D3-brane etc. These ISD 
fluxes stabilize the position of D7-branes and D3-branes [29], when the back 
reaction from D3-branes to the background is neglected [30]. The other open 
string moduli of D7-branes cannot be stabilized in the ISD flux background. 
SUSY breaking effects (including imaginary anti self-dual fluxes) or other 
bulk geometry may stabilize them. 



2.2 KKLT model 

In this subsection, we review on the KKLT model. See Appendix A for 
details of potential analysis. We assume that the dilaton and complex 
structure moduli are frozen around Mstring through flux compactification 
as mentioned in the previous subsection. Here, for simplicity wc neglect 
constant Kahler potential of frozen dilaton and complex structure moduli 
because those only change overall magnitude of scalar potential. To stabilize 
the remaining Kahler modulus, the T-dependent superpotential is added in 
Ref. [4]. Such superpotential can be generated by non-perturbativc effects 
at a low energy scale in the ISD flux background with G(o,3) 7^ 0- Then 
Kahler potential and superpotential are written as 

K = -31n(r + T), (15) 
W = u;o-Ce-"^, (16) 

where a = Stt'^/N with N eN and , 

wo = { Wfiu^ ) = ( <^(o,3) A ) . (17) 



6 



The second term in R.H.S. of (16) originates from SU (N) (N > 1) gauge 
group gluino condensation on non-magnetized D7-branes wrapping unwarped 
4-cycle or Euclidean D3-brane instanton (N = 1) wrapping a similar cycle. 
The constant C can depend on complex structure moduli, i.e. C = C{{U")). 
For gluino condensation, we have C = NM^^^^^g and for D3-brane instan- 
ton C is 1-loop determinant of D3-brane mode which depend on complex 
structure moduli. Thus, the natural order of C is of C(l) or suppressed by 
one-loop factor. This superpotential generates the following scalar potential. 



Vf 



-,K 



\DtW\ 



2{T + Tf 



3|VF|^ 



(18) 



and stabilizes the Kahler modulus, such that StVp = DtW = 0, i.e. 

-3 



DtW 



Wq 



(T + T) 

That corresponds to the SUSY AdS vacuum 

(F^) = -3(e^|T^P). 
Then, we require \wo\ <C |C| and the VEV of T is obtained as 



0. 



aReT = In 



|C|(3 + Q(T + r)) 



ln(?i;o) I — In 



almT = -Arg{wo) + Arg{C) + 2Trn, n G Z, 



(19) 

(20) 

(21) 
(22) 



that is, aReT = 0(47r^). The gravitino mass 1713/2 and the modulus mass 
rriT are obtained as 



m 



3/2 



= {e^\W{') 



~ \wo\ < 1, 



(23) 



mT^a{T + T)m3/2- (24) 

In order to obtain suppressed gravitino mass, we need a very small value of 
Wq, that is, (|W/i„^|) <^ 1. Study on landscape of flux vacua suggests that 
the number of the vacua N^ac with gravitino mass ^71,3/2 is typically given by 
^vac ~ ^3/2 ■lO^'^^ "^2]. This number is amazingly large, and we can tune 
fluxes such that |u;o| ^ 1- For example [6], it is tuned as |G(o,3)/G(2,i)| <^ 1- 
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To realize a dS or Minkowski vacuum, we need to uplift the potential 
by 3777,3^2 ~ l^oP ^ 1- That has been done in Ref. [4] by adding a single 
D3-brane stabilized on the tip of warped throat as an origin of uplifting 
scalar potential. The following potential 

_ y^T^ 1 ^ D 

~ gi 47r2 (i?eT)2 " (T + Tf ' ^ ' 

where = (27r)~^(a')^^ is D3-brane tension, is generated from D3 ten- 
sion and Wess-Zumino term. Because of the warp factor ao, the uplifting 
potential Vl can be very small. Then the total scalar potential is written as 

V = Vp + Vl. (26) 

In order to have {V) = 0, wc have to tune the warp factor using fluxes, such 
as ao = exp[-27r/i/(35'5/)] ~ y/m^ ■ 

The uplifting potential shifts the minimum. The VEV of T slightly 
changes and non-vanishing is generated. Prom the above superpotential, 
we can estimate Wtt = —aWx- Then, by use of analysis in Appendix A, 
we evaluate as 

(27) 



(T + T) 



aReT 



Arg{{F^)) = Arg{{W)) = -Arg{wo). (28) 
It is useful to use the following parameter a. 



a = 



F^ln(Mp/m3/2)' ^ ^ 



in order to represent the ratio between anomaly mediation and modulus 
mediation. In this case we have a ~ 1. We can delete phase of superpotential 

due to U(\)ii symmetry W e^'^^'^'9{wQ)^i and PQ symmetry aT aT + 
i{Arg{C) — Arg{wQ) ) [6, 33]. The results of (23) and (24) do not change. 

We show an illustrating example in Pigure 1. We take parameters as 
C = N, a = Stt^/N, = 5 and L» = 6.3 X lO'^'^. Then, we obtain 
(ReT) ~ 2.2 and 7713/2 ~ 25 TeV. A height of bump at ReT ~ 2.4 in this 
example is 3m3^2 5 because wc uplifted the potential by 3rn^/2- general, 
the height of bump is estimated as 0{my2)- 
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Figure 1: A potential plot of KKLT model with parameters, wq = 10~^^, 
C = N, a = Stt^/N, at = 5 and £» = 6.3 X IQ-^^. Then, we obtain 
{ReT) ~ 2.2 and 171^/2 — 25 TeV. A horizontal axis is t = ReT and a 
vertical axis is V. 



3 Moduli mixing in gauge coupling 

In several string models, the gauge kinetic function is obtained as a linear 
combination of two or more moduli fields. In this section we show concrete 
examples. 

First, in heterotic (M-)theory one-loop gauge coupling is given by [18] 



f strong — S — f5T + /ms, 



(30) 



where /strong is the gauge kinetic function of a strong gauge group and J is 
the Kahler form on CY with hi^i = 1. This can be seen from lOD Green- 
Schwarz term /jv^^p B2 A Ag or IID Chern-Simmons term /^^^^ C3 A G4 A G4. 
The third term /m5 denotes a contribution from M5-brane position moduli 
Y in the orbifold interval, such as /m5 = aY'^/T, a ~ Jqy J /\*&J- Hence, 
gluino condensation may generate moduli mixing superpotential. 

Wgc - ^M-afstrongl a = SlT^N for SU{N). (31) 

That is, the fact that moduli mix in the gauge couplings, implies that moduli 
may mix in non-perturbative superpotential. ^ 

^Moreover, M2-brane stretched between M5-brane and strong coupling M9-brane can 
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Type II models such as intersecting D-brane models or magnetized D- 
brane models have gauge couplings similar to those in heterotic models [21]. 
For example, in super symmetric type IIB magnetized D-brane models on 
T^/{Z2 X Z2) orientifold with /i^"/*^ = 1 ^) gauge couplings are given as 
follows, 

fmD7 = \m7\S + \W7\T, 

fmD9 = mgS-wgT for 03/07 System, (32) 

where rUp, Wp (p = 7, 9) G Z are Abelian magnetic flux contribution F 
from world volume and Wess-Zumino term. In this case, Abelian gauge 
magnetic flux F is quantized on a compact 2-cyclc C2 as J^^^ F G Z. Then, 
rUp {p = 7, 9) is given by mj = J^j^^ F AF and mg = /^^g F A F A F. 
On the other hand, wg is the winding number on a wrapping 4-cycle and 
magnetic flux contribution, wg = J^^g *&Jbuik A F up to a mimcrical factor. 
Moreover, is the winding number of D7-brane on the 4-cyclc. Thus, 
signs of mg and wg can depend on magnetic fluxes and SUSY conditions. 
For example in [35], one can find negative mg and wg. In addition T-duality 
action can exchange the winding number for the magnetic number, but the 
result is similar, that is 

/^09 = WgS-MgT for 09/05 system, (33) 

where Wg, Mg G Z, Wg is the winding number on the 6-cycle and Mg is 
the winding number on the 2-cycle and magnetic flux contributions, Mg = 
ImD9 Jbulk A F A F. Here we have neglected numerical factors again. 
The gauge coupling on these magnetized brane is written by 

= \mgReS - wgReT\ . (34) 

9mD9 

generates [34] 

W^M2 ~exp[-(/3T-y)]. 

This is also moduli mixing superpotential. 

® Actually, /{Z2 x Z2) orbifold, whose orbifold twists are 

uj : [z , z ,z ) ->■ [z , - z , - z ), 

have three Kahler forms in the bulk, that is = 3. We identify indices of those cycles 
for simplicity. 
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Note that magnetic fluxes can contribute to RR tadpole condition of 4-form 
potential and 8-form potential [35] [36] and they should satisfy the tadpole 
cancellation condition, 

^ + E E + Qt"'' = 0' 

p=7,9 a=stacks 

5^ N^w; + Qf^^^' = 0, (35) 

p=7,9 o=stacks 

where Q"*'*^^* {p = 3, 7) are contributions of non-magnetized Dp-brane and 
Op-plane and Np is the stack number of magnetized Dp-brancs. In this 
paper, we treat Wp, nip as free parameters because we only concentrate on 
a stack of magnetized D-branes. 

In Type II A intersecting D-brane models, which are T-duals of the above 
IIB string models, the above expressions of Kahlcr moduli change for com- 
plex structure moduli. However, since there are 3-form and even-form fluxes, 
all geometric moduli can be frozen in this type IIA model at low energy as 
a supersymmetric AdS vacuum. 

In orbifold string theory, moduli in twisted sector, the so-called twisted 
moduli M, can exist [37]. These modes can contribute to gauge kinetic 
function on D-brane near orbifold fixed point, 

fop = {S or T) + aM, (36) 

where a is C?(0.1) — 0(1) of parameter, depending on gauge and orbifold 
group. These twisted moduli may stabilize easily due to their Kahler poten- 
tial [38], but make little contribution to gauge coupling because the mod- 
uli are related to orbifold collapsed cycle. Then, we may naturally have 
(M) < 1. 



4 KKLT models with moduli-mixing superpoten- 
tial 

As seen in the previous section, gauge kinetic functions are linear combina- 
tions of two or more moduli fields in several string models. Thus, in this 
section we study the model that the gauge kinetic function for gluino con- 
densation superpotential is a linear combination of two moduli, say S and T, 
but one of them, S, is frozen around Mstring-, e.g. by flux compactification.^° 

^°In Ref. [22], we will study the model that both of S and T remain light. 
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For concreteness, we use the terminology of models based on 03/07 type 
IIB string theory, in particular the gauge kinetic functions (32). However, 
if we can realize the above situation in other string models, the following 
discussions are applicable for such string models. 

After dilaton and complex structure moduli are frozen out in the ISD 
flux background around M strings we consider the following Kahler potential 
and superpotential at low energy, 

K = -ns ln((5) + (S)) - ht ln(r + T), ng, nr G N, (37) 
W{T) = u;(T)±Be^^^, 6 > 0. (38) 

For the first term w{T) in R.H.S. of (38), we study the following three cases 

w{T) = {Jo^An), (39) 

w{T) = e (40) 

or 

w(T) = Ae~^. (41) 

In the second case, the exponential term of (S) can be generated by gluino 
condensation on a stack of Na D3-branes which is far from warped throat 
At any rate, the first and second cases correspond to w{T) = wq = constant. 
Thus, their potential analysis is the almost same. The third case can be 
realized on a stack of Na non-magnctizcd D7-brancs or Euclidean D3-brane 
which is far from the throat. Alternatively, this term can be generated by 
magnetized D7-branes, and in this case we have A = e~^^ Next, we 

explain the second term Be^^'^ in (38). Since gauge coupling on magnetized 
D-brane is given by /{, = m,ijS ib ujfjT (m^, u'6 G N) like (32), strong 
dynamics on those branes may generate a term like e~'^^ . Thus, the term 
e~^^ [b > 0) means that the superpotential can be generated by gluino 
condensation on (non-) magnetized D7-brane by (32). On the other hand, 
the term e^"^ can imply gluino condensation on magnetized D9-brane and 
we may use this potential so far as a gauge coupling is weak, i.e. 

rribReS > WbReT. (42) 

'^^If we have another inflationary waxped throat, such as Oq ~ 10~^, gluino condensation 
on a stack of N D7-branes or D3-brane on the tip of that throat can also generate a (S) 

dependent superpotential wo ~ (ao) e . 
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In this case from (5) and (42), we need the condition ^ > > 1. 

Now, since the dilaton S is stabihzed at very high energy because of 3-form 
fluxes, then we can write / = rribS it w^T — > mb{S) ± wjjT. Thus, at low 
energy, gluino condensation on these magnetized D-branes can generate a 
very suppressed value of \B\, which is given by 

Q 2 

B = Ce-"'b^^^^ , c=^ ior SUiNb), (43) 

i.e., \B\ <^ 1. Wc assume that C is of 0(1) and one may find b = cwb- 
Then, at the final stage we add the following uplifting potential, 

and the total potential is written as V = Vp + Vl. We tune it such that 
V = Vp + Vl = 0. This uplifting potential for Hp = 2 may be induced by 
adding a few D3-branes at the tip of the throat. If there is a magnetized 
D9-brane, presence of a few D3-branes may give non-trivial effects to the 
D9-brane. In such case, we assume that the compact space is orbifolded, 
D3-branes are located at the orbifold fixed points in the bulk and D3-branes 
are located at the fixed point on the tip of a warped throat like [39]. Then 
the position of D3-branes are fixed and mass of open string tachyon may be 
nearly a TeV scale. Hence it may not affect to our model at this scale. In this 
case, we can have exotic matter on D3-brane at a TeV scale. Furthermore, 
we neglect twisted moduli contribution to gauge coupling as (M) <C 1. We 
also assume that we have a larger amount of total number of D3-branes and 
D7-branes in the bulk than D9-branes, in order not to change the geometry 
of [1]. 

Following [6], we consider arbitrary integer of Up to study generic case. 
At any rate, our models are well-defined as 4D effective supergravity mod- 
els. In what follows, we study four types of models corresponding to all of 
possibilities mentioned above. 



4.1 Model 1 

We study the following superpotential, 

W = wo-Be-^'^. (45) 

In this model, the second term is generated by gluino condensation on a 
magnetized D7-brane, where a gauge kinetic function is / = mb{S) + WbT. 
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If \B\ 3> \wo\, then this is similar to the previous KKLT model. Thus, the 
modulus T is stabilized at, 



bReT ~ In 
Therefore, we obtain 



\B\ 
\wo\ 



(S)-ln|w;o| 



In 



m 



3/2 



(46) 



(47) 



where we have used B = Ce-^''''"^>'{S)/Nt ^ c = N^, b = Sir^Wb/Nb, and the 
gravitino mass is obtained as 777.3/2 — wo- 

When we add the uplifting potential (44) and tune it such that V = 
Vf + Vl = 0, then SUSY is broken and F'^ is induced as 



(T + T) 



=7^777: 



"73/2 



77r b{T + T) 



(48) 



Thus, compared with the results in subsection 2.2, becomes larger by 
the factor In = n^^^"^ ' ^^"^ modulus mass becomes lighter by the 



\B\ 

factor In]^, i.e., uit 



b{T + T) 7773/2- Thus, we obtain 

277T 



a ~ 



377„ 



mi,{ReS) 
wt{ReT) 



+ 1 



(49) 



With 77t = 3, 77p = 2, we have < a < 1. Hence, the modulus mediation 
and anomaly mediation are still comparable except the case with \a\ <C 1, 
where the modulus mediation is dominant. 

When \B\{nT + b{ {T + T) )) < 772'|7t;o|, this analysis may not be reliable 
in perturbative description, because of (ReT) < 0. We may need quantum 
or a' correction to Kahler potential [40]. However, for simplicity, we do not 
consider such case here. 

When the constant term wq is generated by fluxes, we have to fine-tune 
fluxes as G(o,3)/G(2^i) ~ 10^^^, in order to realize soft masses at the weak 
scale. On the other hand, when the constant wq is generated as wq = 
g-87r {S)/Na^ j^Q^ need such fine-tuning for (S) = 0{1). (See also 

Ref. [41].) 

Generic form of potential in this model is similar to Figure 1. The height 
of the bump is of 0(7773/2), and this potential has the runaway behavior at 
the right of the bump. 



14 



4.2 Model 2 



We study the following superpotential, 

W = Ae-"-^ - Be-'''^ . (50) 

This superpotential can be generated by gluino condensation on magnetized 
D7-brane and non-magnetized D7-brane. We assume that (^(0,3) = 0, then 
Wo = for simplicity. However, the following results can change quantita- 
tively but not qualitatively for the case with wq ^ [6] . This model is called 
the racetrack model [42]. For aReT, bReT = 0{4tt'^), the SUSY vacuum 
leads to 



\A\a 

In 

a — b 

nibNa 



ReT ~ In 



B\b\ 

{ReS), (52) 



Nb - WbNa 

"""^^ VmS)^,,J^^n, (53) 



Nb - WbNa ' ' 2{Nb - WbNa) 

ab{T + Tf 

rriT ^ m3/2. (54) 

riT 

Here, for simplicity, we have assumed that ReS ~ 1, a = 87:'^ /Na, b = 
Sir^Wb/Nb, A = Na, B = Nbe-'^>'{S)/wb = ^(i) p^.^^ (5) ^an 

find the physical scales 



2^{ReS) 



Mstring '\Nb-WbNaJ MrK \Nb-WbNa, 

We require the condition that the compactification scale is smaller than the 
string scale, i.e., 

WbNa <Nb< {nib + Wb)Na. (55) 

Then, we obtain 

Wtt , (T + T) 

— — ~ -ab , 56 

Wt riT 

{Ref) = uibiReS) + Wb{ReT) ~ "^"^^ ^r^s). (57) 

Nb - WbNa 
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The gravitino mass is estimated as 



m3/2 c± e 



-a{ReT) 



(58) 



When we add the uplifting potential (44) and tune it such that V = 
Vp + Vl = 0, then SUSY is broken. Generic form of potential is similar to 
Figure 1. By generic analysis in Appendix A, we obtain 



(T + T) 



m 



3/2 



ab{T + r)2 



Since we have 



In 



~ a(ReT) = —(Ref)= ^^'^^^ 



_"i3/2 

we estimate the ratio parameter, 

2 6(T + T) 



(59) 



{ReS), (60) 



a 



rir, 



0(47r2). 



(61) 



That implies that anomaly mediation is dominant in SUSY breaking. These 
results have already been pointed out in [6], but in our model we do not 
need to tune fluxes in order to obtain ReT ~ 1 so far as satisfying (55) and 
ReS ~ 1. Furthermore, since in this case we have 



\W\ ~ \C\ exp 



S-K mil 

Nb-WbNa 



{ReS) 



(62) 



a very small scale may be generated by a value of ReS. That is different 
from the usual racetrack model. 

When we consider the case that both terms in the superpotential are 
generated by magnetized D7-branes, and replace as A Nae~^'^ ma(S)/Na 
with a = Sn'^Wa/Na, then we obtain 



(ReT) 



{Ngmh - Nbmg] 
{WaNb - WbNa) 



(ReS), 



(63) 



and 



{Refa) = ma{ReS)+Wa{ReT)=Na.^;^^^^^^^^^{ReS), 



{WaNb - WbNa) 



(64) 
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In this case, the gravitino mass is estimated as 
and the parameter a is estimated as 

167r2 



a ~ 



WbiT + T) 



o / ma{ReS) , -. \ 
•^^P \wa(ReT) + 



(65) 



(66) 



Its natural order is of 0{4Tr^). That is, the anomaly mediation is dominant. 
In a special parameter region, we may have a = 0{1), where the modulus 
mediation and anomaly mediation are comparable. 

4.3 Model 3 

We study the following superpotential. 



W = Wo + Be 



bT 



(67) 



Now we consider the case that the gauge kinetic function is written by 
fb = i^biS) — WbT, and the condition of (42) should be satisfied. In this 
case, we assume the presence of gluino condensation on a magnetized D9- 
brane. For bReT 3> 1, the modulus T is stabilized at 



bReT 



In 



\wo\nT 



\B\ (b{T + T)-nT)_ 
bImT = Arg{wo) - Arg{B), 



~ In 



\wo\ 



\B\ 



(68) 
(69) 



where wc have used b = 87:^Wb/Nb, B = Ce-^'^»^^^/'"\ C = Nb- One might 
think that for \wo\ S> |-B|, we do not need to tune fluxes to realize ReT ~ 1 
when mb ^ Wb and ReS ^ 1. However, in order to obtain a weak coupling on 
the magnetized D-brane at the cut-off scale, we need the following condition. 



(Refb) = mb{ReS) - WbReT ~ 



Nb_ 
'87r2 



In I Wo I 



1. 



(70) 



Then, we must tune the parameter as \wq\ ~ exp[— Svr^/iV;,]. In the case with 
wq = {J G3 A O ), we have to fine-tune fluxes. However for wo = e~^'^ {S)/Na 



with G 



(0,3) 



0, we do not need such fine-tuning. In this model, since moduli 
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can be stabilized as wq ~ {Be^^), the gravitino mass is given by 7713/2 ~ \u)o\- 
For rub > 1, we can find 



Stt^ 



{Refb) In 



Mr, 



m 



3/2 



(71) 



from (70). 

When we add the uphfting potential (44) and tune it such that V = 
Vp + Vl = 0, SUSY is broken and is induced. Since Wtt = +bWT we 
obtain 



3 



(T + T) '^nTb{T + T) 

2nTWbReT 2nT 



a ~ — 



3np{Refb) 



q / mi{ReS) 
\ wb{ReT) 



1 



(72) 

(73) 



That is, the ratio parameter a is negative. Naturally we would have \a.\ = 
0{l). 

This potential does not seem to have the runway behavior. That has 
important implications on cosmology as will discussed in the next section. 
These properties are remarkably different from the KKLT model, although 
this model only changes to Wb — ^ —Wb from model 1. 



4.4 Model 4 

We study the following superpotential, 

W = Ae-"'^ + Be^^. (74) 

In this model, we assume that non-pcrturbativc effects on the magnetized 
D9-brane and the non-magnetized D7-brane generate the above superpoten- 
tial. We also assume vanishing 3-form flux G(o,3) = 0. In this model, the 
condition (42) should be satisfied. This superpotential can be obtained also 
in heterotic M-theory [43] , where the first term can be originated from mem- 
brane instanton and the second term is originated from gluino condensation 
on the strong coupled fixed plane^^. In this heterotic model the orbifold 
interval T can be stabilized, because of signs of exponents. 

^^See also [44]. 
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This type IIB model is the same as the model 2 except using magnetized 
D9-brane. For aReT, bReT » 1, the modulus T is stabilized as 



ReT ~ In 

a + b 



\A\a 



\B\b 



^""^^ -{ReS), (75) 



Nb + WbNa 



^ ArgiA)-ArgiB) ^_ N^rub 

a + b Nb + WbNa^ " ^ ' 

where we have assumed the same parameters as model 2. One may find that 
ReT tends to be smaller than model 2, and we may need large magnetic 
fluxes m,b- Prom (5) we can find the physical scales 

-2^^{RcS)( ^"""^ Mstrin^_ f NaUlb 



M string '\Nb + WbNaJ MkK \Nb + WbNa, 

Thus, we require the following condition like model 2, 

< iVft < Nairub - Wb). (77) 
We have the gauge kinetic function on the magnetized D9-brane 

{Ref) = "^fe^?' ^R^s). (78) 

Hence, the gauge coupling on the magnetized D9-branc can be weak so far 
as satisfying (77) and {ReS) ~ 1. Since the dynamical scale on this D-brane 
is obtained as |C| exp[— 87r^mb(-ReS')/(7Vb + WbNa)] <^ 1, we can generate 
a small scale as model 2. Furthermore, without fine-tuning, we can have 
moderate values of {ReT) and the gauge coupling so far as satisfying (77). 
Other properties are obtained from model 2 by replacing Wb —Wb- For 
example we have negative a. 

When we replace A ]y^g-s-^^ma{S)/Na ^j^j^ ^ ^ Sir^Wa/Na, results are 
also obtained from model 2 by replacing Wb — ^ —Wb- For example, we obtain 
the parameter a. 



Wb{T + T) 



\ wa{ReT) 

It is negative and its natural order is 0(47r^). 

Here we show an example of potential of model 4 in Figure 2. It is quite 
different from model 1 and model 2, while in model 3 a similar potential 
is obtained. That has cosmologically important implications, as will be 
discussed in the next section. 
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ln[V] 




t 



Figure 2: The potential of model 4 with parameters, ut = 3, ns = 1, np = 
2, = 4, Nh = 5, ruh = 6, Wb = 3, {ReS) = 1, D = 2.2x lO'^^. Then, we 
obtain (i^eT) ~ 1.4, M string /Mp ~ 0.22 and mg/s 40 TeV. A horizontal 
axis ist = ReT and a vertical axis is ln[V]. This potential can make sense 
until ReT < mi,{ReS)/wb = 2. 



5 Implications on SUSY phenomenology and cos- 
mology 

In this section, we discuss implications of our results on SUSY phenomenol- 
ogy and cosmology. 

Concerned about SUSY breaking, our results are qualitatively similar 
to those in Ref. [6], that is, in model 1 and 3, the modulus mediation and 
anomaly mediation are comparable, while in model 2 and 4, the anomaly 
mediation is dominant. Our models generalize those results. Since the 
parameters A and B can expand wider range in our models, the parameter 
a varies a wider parameter region, even if ut and Up are fixed. In model 3, 
we obtain a negative value of a, which would be naturally of 0{1). In model 
1, we have < a < 1 when we fix n^- = 3 and Up = 2. In addition, a value of 
a in model 2 and 4 is of 0(47r^). However, whether the modulus mediation 
and anomaly mediation are additive or destructive depends on gauge kinetic 
functions of visible sector. Suppose that the gauge kinetic function of the 
visible sector is given as 

= m^S + WyT, (80) 
where Wy can be positive and negative. Then, whether the modulus media- 
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tion and the anomaly mediation are additive or destructive depends on the 
signs of a and w^. Furthermore, the size of gaugino masses induced by 
is written as 



1/2 



^ 1 



(81) 



Thus, it can be enhanced or suppressed. At any rate, soft SUSY breaking 
terms in our models have a rich structure. We will study their spectra and 
phenomenological aspects elsewhere. 

Next, we discuss cosmological aspects of our models. The anomaly me- 
diation is sizable except the special case with |a| ^ 1. That implies that 
soft masses in the visible sector are suppressed by one-loop factor compared 
with the gravitino mass m3/2, that is, the gravitino can be heavier like 
^3/2 — C'(IO) TeV. The modulus can have much larger mass. Therefore, we 
may avoid the gravitino problem and the moduli problem in all of models 
like the usual KKLT model [6, 9, 10, 11]. Furthermore, the potential forms 
of model 3 and 4 have more interesting aspects as discussed below. 

The potential of model 1 and 2 has the runaway behavior at the right 
region of the bump like Figure 1, and the height of the bump is of 0(771,3^2)- 
When an initial value of T is in the right of the bump, T goes to infinity. 
Furthermore, when an initial value of T is quite small T ^ 1, T overshoots 
the favorable minimum and goes to infinity. Thus, we have to fine-tune the 
initial condition such that T is trapped at the favorable minimum. That is 
the overshooting problem [45]. 

This type of potential has another problem, that is, destabilization due 
to finite temperature effect. The finite temperature effect induces the addi- 
tional potential term [46], 



where T denotes the temperature. The coefficients, ao and 02 , are written 
by group factors of massless modes, and in most of case 02 is positive. For 
example, we have 03 = ^{^l ~ l)(A^c + 3iV/) for SUSY SU{N^ gauge 
theory with Nj flavors of matter multiplets. In models 1 and 2, the gauge 
kinetic function is obtained as / = mS -\- wT. Thus, the potential term due 
to finite temperature effect is written as 



Ay = (ao + a29^)f\ 



(82) 



AV = ao + a2 



1 



(83) 



m{S) + wT 
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This term destabilizes T at not so high temperature, but the temperature 
corresponding to the intermediate scale [47] . 

The above two problems arc not problems only for our models 1 and 2, 
but are rather generic problems. On the other hand, the potential in models 
3 and 4 has the term like Be^-^ . Such term may avoid the overshooting 
problem, and this term is reliable except when Be^^ > 0(1), because in 
that region we would have uncontrollable effects. However, this reliable 
region of the potential is much higher than 0(171^2)1 which is the height 
of the bump in model 1 and 2 as well as other potentials with this type of 
potential forms. The same behavior of the potential is helpful to avoid the 
destabilization problem due to finite temperature effects. For example, in 
model 3 the gauge kinetic function is written as / = mS — wT. Thus, the 
potential term due to finite temperature effects is written as 



AV 



1 



m{S) — wT 



T". (84) 



That makes T shift to a smaller value, because smaller T corresponds to 
weaker coupling. Therefore, the VEV of T does not destabilize. Model 4 
also has the same behavior. As results, the potential form in models 3 and 4 
are cosmologically interesting from the viewpoint to avoid the overshooting 
problem and destabilization problem due to finite temperature. 



6 Conclusion 

We have considered the KKLT model with moduli-mixing superpotential, 
assuming that one of them is frozen. Such superpotential can be obtained 
e.g. by gluino condensation on magnetized D-branes, while it may be gener- 
ated in other setups. We have studied four types of models. In these models, 
the hierarchy between the Planck scale and gravitino mass can be written 
by gauge coupling, such as ln(Mp/m3/2) — 8Tr'^{Ref) /N, that is, magnetic 
fluxes can generate a large hierarchy. Model 1 is almost the same as the 
KKLT model, but the ratio a between anomaly mediation and modulus me- 
diation can take various values. Models 2, 3, 4 do not require fine-tuning of 
3-form fluxes to realize {ReT) 1 because of very small coefficient B. How- 
ever, model 3 needs fine-tuning of 3-form fluxes in order to obtain a weak 
coupling on the magnetized D9-brane, but in the case that wq can be gener- 
ated by gluino condensation on D3-brane, we may not need to tune minutely. 
We may need to tune slightly the open string sector such as magnetic fluxes. 
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the winding number and the number of D-branes instead of 3-form fluxes 
(in the closed string sector). In model 3 and 4, a becomes negative. All of 
models lead to a rich structure of SUSY breaking including new patterns of 
soft SUSY breaking terms. Such spectra and their phenomenological aspects 
would be studied elsewhere. 

In most of models except for |a| <C 1, the gravitino and moduli masses 
are of 0(10) TeV or much heavier. Such spectrum is important to avoid 
the gravitino problem and the moduli problem. Furthermore, the potential 
form of models 3 and 4 have good properties for cosmology because of the 
exponential factor, exp[+6T] for 6 > 0. That may avoid the overshooting 
problem and the destabilization problem due to finite temperature effects. 

We have studied the models that two moduli fields S and T have mixing 
in superpotential, assuming one of them S is frozen already around the string 
scale M string- Alternatively, both of them may remain light. In Ref. [22], 
we will study such models, i.e., moduli stabilization and SUSY breaking. 
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A Analysis on the KKLT type scalar potential 

In this appendix, we summarize analysis on the KKLT type scalar potential. 
A.l SUSY potential 

First we consider the following F-term scalar potential Vp of supergravity 



model. 



Vp = GiG-jG 
G = K + hi\W\^ 



- 3 



(85) 
(86) 



where Gi = djG, Gj-j = djdjG. 
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The first derivative of Vp is obtained as 

+ Gj; ■ [GkiG'''^ + Gk (S/G^^)) 



= F 



K 



+ F 



Kj^-^diF^ + Se^/'^WK 



IL 



(87) 



(88) 



where we have defined F^ = -K^-^ e^/'^DjW . Thus, the SUSY point, i.e. 
Gj = DjW = 0, satisfies the stationary condition, diVp = 0. At such SUSY 
point, mass matrices are given by 



Vij = -e^Gij, 
Vn = -e^Gjj, 



IJ 



2Gjj + GikG-jyG 



KL 



(89) 
(90) 

(91) 



Here we concentrate to the model with only single field X. In this case, 
mass matrix is given as 



V^j^ + Re{Vxx] 
-Im{Vxx) 



-Im{Vxx) 
Vxx-MVxx) 



(92) 



in the basis of {Re{x), Im{x)), where x = X — Xq with Xq satisfying 
dxV\x=Xo = 0- Their eigenvalues are obtained as 



XX \ 



(93) 



Therefore, the SUSY point corresponds to the minimum of the potential if 
^xx > \Vxx\, that is. 



\Gxx\ > "^Gxx 



(94) 



A. 2 KKLT type potential with uplifting potential 

Now we consider the KKLT type potential. We use the following form of 
Kahler potential and generic form of superpotential. 



K = -utHT + T), 
W = W{T). 



riT e N, 



(95) 
(96) 
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The original KKLT model has nr = 3. The SUSY point, i.e., DtW = 0, is 
obtained as 



(T + T) = nr 



W 
Wt 



e R. 



(97) 



This point corresponds to the minimum of Vf if the VEV of T satisfies 



nT{nx — 1) Wtt 



(T + r)2 



w 



> 



2nT 



{T + Tf ■ 



(98) 



We consider the case that the above condition is satisfied, but such minimum 
of Vp is the SUSY AdS vacuum with the vacuum energy Vp = —3e^\W\^ = 
— 3771,3^2- 

In order to realize the dS vacuum, we add the following form of uplifting 
potential, 

Vl = 7—^—, < L>< 1. (99) 

In the original KKLT model, this is originated from D3-bane and we have 
np = 2. Here following Ref. [5, 6], we consider generic integer for np. Now 
the total potential is written as 



(100) 



and adding Vp changes only vacuum value of ReT from (97) but not ImT. 
We demand that a change of ReT is small and the cosmological constant 
{V) vanishes, that is. 



{VL) = -{VF)c:,3ml/^. 

Prom (88), the first derivative of V is written as 



(101) 



(T + r)3 



(T + Tf 



— T 



riT 



(T + r)2 



TT7 riT 



(T + T)2 



{T + T) 



Vl. 
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Now, we expand as T = TgusY + 5T such as DtW\t=Tsusy — 0- Then, the 
F-term and its derivatives are evaluated as 

pT ^ {T + T) _ t)Wt), (102) 

riT 

drP'^ ~ a{T + T)-e^/'^W = a{T + T)m^l2, (103) 
~ -e^/'^W = -7113/2. (104) 
Here we have used Wt — nrW/ (T + T) and defined a real parameter a, 

We have also assumed that |a| • ReTguSY S> 1. There is an important point 
for CP phase. Since 

ST, ^, 5^eR. (106) 

the CP phase of W and F arc the same. That is, when wc write W = 
\W\e^^^ , then we have F"^ oc e'*^*^. Using the above results, we can write 

OtV ^ m3/2^^T^-[a(r + T)F^]-^^;^l^^. (107) 

Therefore, for aReT ^ 1 the stationary condition drV = and the condi- 
tion for the vanishing cosmological constant lead to 

- n,^^^e-^^^, (108) 

- ^ < 1. (109) 



TsuSY 2 nra'^iReTsusY? 

Furthermore, the F-component of conformal compensator superfield (f) is 
given by 

E± = e'^'^W + \kiF' , (110) 

90 3 



e-''^^m3/2, (111) 

where (pQ (g R) is a scalar component of (j). Hence, CP phases of and 
F'^ are aligned [6, 33]. It is useful to define the ratio a as follows, 

F^{T + T) _ nTa{T + f) 
00 ln(Mp/m3/2)i^^ 3npln(Mp/m3/2) ' 

because of ln(Mp/m3/2) ~ C'(47r^). 
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